Gas density profile in dark matter halo in chameleon cosmology 
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We investigate the gas density, temperature, and pressure profiles in a dark matter halo under the 
influence of the chameleon force. We solve the hydrostatic equilibrium equation for the gas coupled 
with the chameleon field in an analytic manner, using an approximate solution for the chameleon field 
equation with the source term, with a generalized Navarro-Frenk-White universal density profile. 
We find that the gas distribution becomes compact because a larger pressure gradient is necessary 
due to the additional chameleon force. By confronting the theoretical prediction with the data of the 
temperature profile of the Hydra A cluster according to Suzaku x-ray observations out to the virial 
radius, we demonstrate that a useful constraint on a model parameter can be obtained depending 
on the value of the coupling constant. For example, the upper bound of the background value of 
chameleon field, (j)oa < 10~''Mpi, is obtained in the case (3—1, where /3 is the coupling constant 
between the chameleon field and the matter, and Mpi is the Planck mass. However, the error of 
the present data is not so small that we obtain a useful constraint in the case /3 = l/\/6, which 
corresponds to an f{R) model. 
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The accelerated expansion of the Universe is one of 
the most fundamental mysteries in basic science. Long- 
distance modification of the gravity theory is a challeng- 
ing approach to this problem. However, any gravity the- 
ory must pass the stringent constraints from the Solar 
System. The chameleon mechanism is a noble mecha- 
nism for screening a scalar degree of freedom which ap- 
pears in a class of modified gravity models, depending 
on the density of matter in the local environment [ij, y] . 
Newtonian gravity is recovered in a high-density region, 
thereby evading the Solar System constraints. Recently, 
it has been pointed out that the modification of grav- 
ity might be detected using halos of galaxies and galaxy 
clusters, because the screening mechanism could not be 
complete in their outer regions |3H8j . 

In the present paper, we focus on the gas distribu- 
tion in a dark matter halo under the influence of the 
chameleon force. In Ref. [9j, the authors found an ana- 
lytic solution of the chameleon field, assuming the matter 
distribution of the Navarro-Frenk-White (NFW) univer- 
sal density profile [lOJ (c.f. Ref. [11|). Utilizing their 
analytic method, we investigate the gas density, temper- 
ature, and pressure profiles under the influence of the 
chameleon force. We find that the chameleon force sig- 
nificantly influences the gas distribution. We also demon- 
strate a useful constraint on the chameleon gravity model 
from confronting the theoretical temperature profile with 
x-ray observations of a cluster of galaxies. 

The chameleon field equation for a quasistatic system 
in the Einstein frame is 
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where V is the potential, p is the matter density, /3 is 
the coupling constant, and we have defined the reduced 



constant G. Here, we assume F(0) = A'*+"/0", where 
A is the mass dimension parameter and n is the dimen- 
sionless parameter. We also assume I3(p/Mp\ ^ 1. The 
coupling between the scalar field and the matter density 
is the key for the chameleon mechanism, as we see below. 

We follow the analytic method in Ref. [9| to find a 
solution for Eq. ([T|). In the present paper, we assume 
the generalized NFW density profile p{x) = ps/x{l + x)'' 
with X = r/vs, where p,, and r^ are the character- 
istic density and scale of a halo, respectively, and b 
is a parameter. The NFW density profile is the case 
b — 2. The mass within the radius x of the halo is 
given by AI{x) = Airr^ J^ dxx^p{x). Instead of the pa- 
rameters Ts and ps, we introduce the virial mass Afvir 
and the concentration parameter c, which are defined by 
A/vir = (47r/3)rvir'^AcPc and c = r^ir/rs, where Ac is 
the ratio of the spherical overdcnsity p{< c) within the 
virial radius r^ir to the critical density of the universe Pc] 
i.e., Ac = p{< c)/pc, for which we adopt Ac = 100 in a 
spherical collapse model. 

The analytic solution for Eq. ([T|) is obtained by match- 
ing the interior solution <j)\nt and the exterior solution 
4>ont, where (/)int is given by solving Ep. ([1]) while neglect- 
ing the term on the left-hand-side, while 4>oMt is given by 
neglecting the first term on the right-hand-side of Eq. ([T]). 
Then, we find 
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where we have defined 0s = ('T-A"+^Mpi//3/9s) 



i/("+i) 



Planck mass by Mpi = l/(87rG') with the gravitational B = (Spsr'^/Mpi, and C and Xc are determined by solving 



the matching conditions at x 



equilibrium gives 
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The validity of the analytic solution is demonstrated for 
the case 6 = 2 in Refs. ja, l9|. 

Note that (ps is the typical value of the chameleon 
field in the interior region, where the chameleon mech- 
anism works. This means that 0s ^ 0oo, because the 
chameleon mechanism screens the chameleon field. With 
this fact, Eqs. (|3]) and Q are approximated as 



C-B{{l + x,) 



2-6 



\l-6 



l)/ib^2) + (f>^xc, (5) 
, 0, (6) 



since we are considering the case Xc = 0{1). Hence, the 
scalar field in the exterior region is independent of n and 
A. This is important because the constraint we obtain 
becomes independent of n and A. 

In our modeling of the gas distribution in a dark matter 
halo, we make a few assumptions for simplicity. First, we 
assume that the dark matter dominates the dark halo de- 
scribed by the generalized NFW density profile. The halo 
density profile could be affected by the modification of 
the gravity theory; however, we assume the same profile, 
since its validity is partially supported by iV-body simu- 
lations for the Dvali-Gabadadze-Porrati model and f{R) 
model [SHsl, |l2| . Recently, Ref. [5| confirmed this validity 
for an f{R) model, and argued a qualitative explanation 
for the validity. Second, we assume that the baryon den- 
sity is negligible in the dark matter halo, which allows 
us to neglect its effect on the scalar field equation. The 
effect of the baryon component is discussed in Ref. [13[ . 
The validity of this assumption is also supported by the 
recent measurements of the density profile of a cluster 
halo through gravitational lensin g, w hich show that the 
NFW profile fits the data well fllJl5l|. Third, we assume 
that the scalar field is coupled with the baryon compo- 
nent as well as the dark matter component. For example, 
in an f{R) model, the chameleon force is coupled with 
both the dark matter and the baryon components. 

Within Newtonian gravity, a useful model of the gas 
density profile is considered in Refs. 13|, |l6|. By assum- 
ing hydrostatic equilibrium between the gas pressure and 
the gravitational force from the dark matter with the 
NFW density profile, the universal gas density, temper- 
ature and pressure profiles are derived. We follow the 
method of Refs. 13|, ll6|, but taking into account the 
chameleon force as well as the gravitational force, we de- 
rive the gas distribution in a halo. Now the hydrostatic 
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where pg and Pg are the gas density and the pressure, 
respectively, and e = unless explicitly stated other- 
wise. Here 4>g denotes the gravitational potential, given 
by solving the gravitational Poisson equation, A^g = 
p/(2Mpi). For the generalized NFW density profile, we 
find the solution 0g(x) = (t)o[l - {1 + xf'^ / {b - 2)x\, 
where we define 0o = —psr1/2Mp\{h — 1). We as- 
sume that the gas obeys the polytropic equation of state 
Pg ex PgTg ex p^ with the polytropic index 7 and the 
gas temperature Tg. Introducing the function yg{x) by 
Pg{x) = Pg{Q)yg{x), Pg{x) = Pg{0)y]{x), and Tg{x) ^ 
Tg{0)y2~^{x), we obtain the solution 
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where we have defined A = —pmp(j>o{'j — 
l)/kTg{0)Mpi{l + e)7, k is the Boltzmann con- 
stant, and iirUp represents the mean molecular mass. 
We determine the parameter 7 by Expression (17) in 
Ref. 17[. Our conclusions are not altered qualitatively 
for the assumption on 7 within 1.1 < 7 < 1.3. 

Figure [1] shows the gas density profiles, comparing the 
case with the chameleon force (solid curves) and the case 
of Newtonian gravity (dashed curves), adopting virial 



masses of M^ 



Mr, 



IO^^Mq, W^Mq, and 4 x IO^^i^iq, 
from top to bottom, respectively. The gas density de- 
creases rapidly in the outer region (see the solid curves), 
where the chameleon force is infiuential. For the large 
mass cluster, the chameleon mechanism works out to 
large radii, because the density of dark matter is high 
enough even outside the halo. On the other hand, for 
the small mass cluster, the chameleon mechanism works 
only at small radii, because the dark matter density is 
high only in the central region. Because the chameleon 
force is an attractive force, a larger pressure gradient is 
necessary for balancing between them. This makes the 
gas distribution compact. This feature is more significant 
for the smaller-mass halo. 
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FIG. 1: Gas density profile as a function of the radius r/vs- 
The solid (dashed) curves are with (without) the chameleon 
force, with virial masses Mvir — 4 x IO^^Mq (green light 



curve), 10 Mq (blue dark curve), and Mv; 
curve), from top to bottom, respectively, 
adopted 13 — 1, n — 1.8 x 10" 
0^ = 1.1 X 10~^Mpi. 
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FIG. 3: The contours of Ax^ on the parameter plane 0oo- 
Afvir. Here we have fixed /3 = 1 and b = 2, but c and Tg{0) 
are varied as fitting parameters. The contour levels of the 
inner dashed curve and the outer solid curve are Ax^ ~ 2.7 
and 6.6, respectively. 
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2: Temperature profiles as a function of the ra- 
r±. The points with error bars show the observa- 



FIG. 

dius 

tion data of the filament direction [1^. The curves show 
our theoretical model. The solid curve adopts the best- 
fit parameters {(p^, A'Uir,c,Tg{0)) = (5.4 x 10"^Mpi,5.1 x 
1O^''M0,5.8, 4.9keV). The dashed curve adopts (j)^ ^ 2 x 
10~^Mpi, while the dotted curve adopts (pea = 1.3 x lO'^Mpi, 
where the other parameters are the same as those of the solid 
curve. Here we have fixed /3 = 1 and b = 2. 



Using this characteristic feature, let us consider a con- 
straint on the chameleon gravity model. To this end, we 
consider x-ray observations of a cluster of galaxies. Be- 
cause of the steep drop of the gas density in the presence 
of the chameleon force, a similar drop in the x-ray sur- 
face brightness may appear in the outer region. In the 
present paper, we compare the x-ray temperature pro- 
file from the data reported from Suzaku observations of 
the Hydra A cluster out to the virial radius [18|. The 
Hydra A cluster is a medium-sized cluster located at a 
distance of 230 Mpc. Two different fields arc observed 



in Ref. pjj. One is the northwest offset from the x-ray 
peak of the cluster, and the other is the northeast off- 
set. The former and latter fields are called the filament 
and void, respectively, because each field continues into 
the filament and void structures. In Fig. [51 the points 
with error bars show the data of the filament direction 
in Ref. ^. 

The curves in Fig. [2] show our theoretical model of the 
x-ray surface brightness temperature, computed with the 
formula 
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where Ac(Tg) is the cooling function, for which we have 
assumed Ac(Tg) (x Tg' (e.g., Ref. [l9|), and r± is the 
radial coordinate perpendicular to the line-of-sight direc- 
tion. The solid curve is the best-fit curve, whose param- 
eters are noted in the caption. The dashed curve and the 
dotted curve adopt ^cx, = 2x 10"^Mpi and 1.3x lO'^Mpi, 
respectively, where the other parameters are the same as 
those for the solid curve. 

The dotted curve, the solid curve, and the dashed curve 
in Fig.[2]represent the characteristic curves which appear 
when we vary (poa from a sufficiently large value to a 
smaller one. First, the dotted curve represents the limit 
of the modified gravity. Namely, for the large value of 
^oo ^ B, Xc becomes negative from Eq. (|6]). This means 
that there appears no interior region in a halo where the 
chameleon mechanism works to recover Newtonian grav- 
ity. Thus, for the case (j)oo > B, we have (/)(x) = (l)out{x) 
for the entire region, and therefore the solution Eq. ([5]) 



should be replaced with 
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On the other hand, the dashed curve is the same as the 
limit of Newtonian gravity. For a small value of 0oo , we 
have a large value of Xc from Eq. ^. This means that 
the chameleon force is influential only at very large radii. 
Note that the interior solution yg{x), Eq. ^ for x < Xc, 
can be approximated by taking the limit of /? — >■ in 
Eq. (fTT|) , because (j)s takes a very small value to screen the 
scalar field where the chameleon mechanism works. In 
summary, the dotted curve and the dashed curve are the 
two opposite limits, and our theoretical curve is restricted 
by these two limits. Note that the limit of the modified 
gravity Eq. pT|) depends on the coupling constant /?. 

Figure [3] shows the contours of Ax^ on the parameter 
plane for ^oo and Mvir, where x^ is simply defined by 
X' = Y.l=l{Tx{r^,^) - T°^'-)y{AT°'''-y-,' where T^'^^- 
and AT"'^^- are the observed data and the error of the 
filament direction, respectively, and Tx(r±,i) is our the- 
oretical model. Here, we have fixed /3 = 1 and 6 = 2, but 
the parameters c and Tg(0) are varied so as to minimize 
X^ within the range 3 < c < 10 and T^/a < Tg{0) < J^a 
with a = 1.1, where Tg* is given by Eq. (19) in Ref. ^. 
When taking Tg{Q) as a completely free parameter, it 
is difficult to obtain a useful constraint from the present 
data due to the degeneracy between Tg{Q) and Mvh-. The 
minimum value of x^ is 1.0, where the number of degrees 
of freedom is 3. The behavior of the contour is explained 
by the fact that the theoretical curve approaches that 
of Newtonian gravity as (poo becomes small and that the 
steep drop becomes significant as (jiao increases. Figure [3] 
gives an upper bound of (p^c < 10~^Mpi at the 2-sigma 
level for the case 6 = 2 and /? = 1. We obtain a similar 
upper bound of 0oo for different values of 6, which are 
summarized in Table HI The upper bound of 4>oo becomes 
larger for smaller b, but we may conclude that the re- 
sults do not significantly depend on b. Table [J includes 
the results with the void direction. The upper bound of 
4>oo depends on the data, i.e., the filament direction and 
the void direction; however, our conclusion does not alter 
significantly. 





Upper limit for (j>aa in unit of [AIpi] 




Filament 


Void 


6=1.7 


1.4 X 10"'' 


0.9 X 10"* 


6 = 2.0 


1.0 X 10"^ 


0.8 X 10"* 


6 = 2.5 


0.8 X lO^'' 


0.6 X 10"* 



TABLE I: Upper bounds of 0oo at the 2-sigma level for dif- 
ferent values of b and the data for the filament and void di- 
rections. Here we have fixed 13 = 1. 



So far, we have considered the case /3 = 1; let us now 
discuss the case /3 = 1/VQ, which corresponds to an f{R) 
model [20|. In this case, we could not obtain a useful 
constraint on (/)oo, which is explained as follows. The 
theoretical density profile is limited by two characteris- 
tic curves, Eq. ^ and Eq. ^ with /? = 0. When /3 
is small, the difference between these two characteristic 
curves is small, because the drop of the gas density is 
not steep. This is the reason why no useful constraint 
on the f{R) model was obtained from the present x-ray 
data here. 

In summary, under the assumption of hydrostatic equi- 
librium between the gas pressure gradient and the gravi- 
tational and chameleon forces, we derived the gas density 
profile in an analytic manner. Here we assumed the poly- 
tropic equation of state for the gas and the generalized 
NFW density profile for the dark matter distribution. 
The chameleon force may give rise to a steep drop in the 
gas distribution in the outer region of a halo. This fea- 
ture is more significant when the mass of a halo is small 
and /? and 0cx3 are large. The gas density profile depends 
on /3 and 4>oo, but it does not depend on n and A. This 
provides us with an opportunity to constrain /3 and (/jqo 
by comparison with observations. We demonstrated a 
constraint on (jjao in the chameleon gravity model, us- 
ing the data of the temperature profile from the x-ray 
observation [18|. We obtained a useful upper bound of 
4>oo < 10~'^Mpi in the case /3 = 1 [2l|; however, no useful 
constraint was obtained in the case /3 = l/v6, which cor- 
responds to an f{R) model. In order to obtain a useful 
constraint, observations of the outer region of a smaller 
mass cluster are more advantageous. Furthermore, a 
combination with other observations like the weak lens- 
ing measurements might improve the constraint. In our 
investigation, the assumption of the hydrostatic equilib- 
rium of hot gas might be crucial. To estimate the effect 
of deviation from it, we obtained similar constraints by 
adopting the nonzero values of e = ±0.5 in Eq. ([7]). The 
upper bound of (f>oo changes from lO^'^Afpi for e = to 
0.6 X lO-^Mpi and 2.1 x lO'^Mpi for e = 0.5 and -0.5, 
respectively. Thus, the assumption of hydrostatic equi- 
librium is crucial to the constraint, but we may obtain 
a useful constraint if we can model the state of the gas 
correctly. Further study is necessary for this problem. 
Finally, we assumed spherical symmetry for a halo, an 
assumption whose validity must be checked when com- 
paring with observational data. In the present paper, the 
results in Table U do not depend on the filament direc- 
tion and the void direction significantly, which suggests 
the validity of our assumption. 
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